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1. Introduction 

The calculation of the relevant thermodynamic potential of an interacting many- 
body system in the thermodynamic limit belongs to the most fundamental problems 
of statistical mechanics. For Yang-Baxter integrable models, very sophisticated 
techniques like the thermodynamic Bethe ansatz (TBA; for a review, see [1]) and 
the quantum transfer matrix (QTM; for a review, see [2]) have been developed which 
yield exact expressions for the thermodynamic potential in terms of certain auxiliary 
functions. These auxiliary functions are determined by non-linear integral equations, 
which are well suited for numerical investigations (see, for example, [3, 4, 5]). 

As far as the spin-1/2 XXZ-Heisenberg chain is concerned, high-temperature 
expansions have been performed analytically (see, for example, [1, 6, 7]) up to high 
orders. At low temperatures, the leading ("universal") terms have been obtained from 
the exact solution. In the absence of a magnetic field, higher orders were calculated 
by combining the T = Bethe ansatz with an effective field theory for the low- 
lying excitations [8, 9]. It still remains a challenging problem to perform such a 
low-temperature expansion directly from the exact solution for the thermodynamical 
potential at finite temperatures. 

On the other hand, for integrable bosonic systems, only the "universal" terms 
stemming from the Luttinger-liquid character are known for general coupling and 
particle density [10, 11]. A virial expansion of the integrable Bose gas was given in [3]. 
Very recently, the low-temperature behaviour of an anyonic gas at large interaction 
strength was considered [12]. However, a systematic low-temperature expansion for 
integrable bosonic systems with arbitrary interaction and occupation has not been 
performed yet. 

In this article, a Sommerfeld-type expansion is applied to the TBA-/QTM- 
equations for integrable systems in an external field, namely the g-Bose model with 
a chemical potential and the critical XXZ-chahi in a magnetic field. The results 
are directly transferable to the integrable Bose gas, and the method itself can be 
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employed for other integrable systems with external fields. The importance of this 
work lies in the close connection between the low-energy effective Hamiltonian and 
the low-temperature expansion of the thermodynamical potential, as pointed out in 
[8, 9] for the zero-field case of the Heisenberg chain. The low-lying excitations of this 
model, as well as of the qSosc model, arc described by a Gaussian model. In [8, 9], 
the leading irrelevant operators and its coefficients were identified for the Heisenberg 
chain in zero field. For finite field, the leading irrelevant operators have been found 
recently [13]. It is an interesting problem now to relate the low-temperature expansion 
found here to the coefficients of the irrelevant operators in this effective Hamiltonian 
for the spin chain. Besides that, for the g-Bose model, only the scaling dimension of 
the leading irrelevant operators are known, but not the operators themselves nor the 
coefficients. This work constitutes a first step to fill this gap. The knowledge of the 
effective Hamiltonian allows for the calculation of asymptotics of correlation functions, 
[8, 9, 13]. 

Our main result is an asymptotic low-temperature expansion of the thermody- 
namic potential g = g{T, n) in even powers of the temperature, 

giT,^i)=Y,g.Ap)T^\ z/ = 0,l,2,... . (1) 

The coefficients gi^in) depend on the chemical potential /i only. For the Heisenberg 
chain, n is replaced by the magnetic field h. We emphasize that Eq. (1) is not a Taylor 
expansion in T, but rather an asymptotic expansion. This means that the coefficient 
in order T^"* is defined as 

The first two coefficients in Eq. (1) are well known: go constitutes the ground 

state energy. Furthermore, gi = — 7r/(6w), where v is the velocity of elementary 
excitations ("universal term"). Our approach allows for the calculation oi gv>2, and 
the calculation of g^ is performed explicitly in this article. 

The universal term originates in the Gaussian part of the effective Hamiltonian, 
with scaling dimension two. The higher-order terms stem from irrelevant operators 
with scaling dimensions 2v [v > 2). This is consistent with the absence of 
backscattering processes in both the g-Bose model and the Heisenberg chain 
in a magnetic field: For the g-Bosons, the conserved particle current excludes 
backscattering [14]. On the other hand, a finite magnetic field applied to the 
Heisenberg chain moves the system away from half filling (in the fermionic picture) 
and thus forbids backscattering at low energies as well. 

This article is organized as follows. In the next section, the low-temperature 
expansion of g for the q-Bosc model is obtained from the corresponding TBA- 
equations. The results are checked for a certain choice of the interaction parameter (the 
phase-model, [15, 16]) by comparing with an independent expansion in this case and 
with numerical results. The third section deals with a high-temperature expansion of 
the q-Boson model, whereas the fourth section contains the low-temperature expansion 
for the critical spin-1/2 XXZ-chain. Calculations not directly necessary for the 
understanding of the main text 'atc d(^f(Trcxi to the appendix. In all calculations, 
units are chosen such that Boltzmann's constant /cb = 1. 
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2. Low-temperature asymptotics of the q'-Bose model 

2.1. Definition of the model 

The g-Bose model is a one-dimensional integrable bosonic lattice model that in 
an appropriate continuum limit leads to the Lieb-Liniger Bose gas. The lattice 
rcgularization stems from a ^-deformation of the underlying commutators [17, 18, 19]. 
Very recently, ground state and thermodynamic properties have been calculated 
numerically and the Gaussian part of the effective Hamiltonian has been identified 
[11, 14]. Let us shortly review the definition of the model. 

g-deformed bosonic annihilation and creation operators are defined as 

bbt = 1 - g2n+2 , btb = 1 - , 

bn=(n-M)b, b^n = (n - l)b^ 
with < g < 1. In the following, it is often convenient to use the parameterization 

g = e-^ r?e]0,oo[. (2) 
A g-oscillator is represented in the Fock space: 

n\n)=n\n); n = 0, 1, 2, 3, . . . ; b|0) = . (3) 

We deal with a chain of length L where the local q'-oscillator algebra is assigned to 
each site i= 1,. . . ,L. The quantities 

L 

l=\ 

1 ^ 1 ^ 

^+ = T3^Eb^bt,+, V-b,bt,_, , 

^ 1=1 ^ e=i 

commute pair-wisely. The eigenvalue of J\f is the number of particles N. The choice 

of the Hamiltonian 

H - -^{V+ + r^)-fiM, (4) 

where is the chemical potential, leads [17, 18, 19] in the continuum limit to the 
Lieb-Liniger Bose gas with Hamiltonian 

n = ^ J [a^*^(x)a^*(x) + 2c'iiHx)¥{x)-i'{x)-^{x) - /x*t(x)i'(x)] dx. (5) 

Here A is the lattice constant of the g-Bose model and the coupling constant c has been 
introduced such that r] = c - A. Furthermore, trivial constants have been absorbed 
into 11. 

Using the TBA, the thermodynamic potential potential (which here is —p, with 
p being the pressure) was found to be given by [19] 

.9(T,M)=-T|^n(l + e-W/^)g, (6) 

where the function e(k) is obtained from 

= - - . - r£ .n (1 . e--) |: . (7, 

These two equations arc valid for arbitrary T. They constitute the starting point for 
our low-temperature analysis. We will calculate contributions to g including O(T^). 
However, higher orders can be calculated as well in the approach sketched below. 
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2.2. Low-temperature expansion 

At T = 0, the function e(A:)|T=o has two zeroes ±B{T = 0) =: ±Bo. In [11], the 
ranges of the particle density n := N / L, the chemical potential ji and the parameter 
-Bo were given as in table 1. Here, we are dealing with temperatures T 4C i.e. 



n 





oo 




-1 





Bo 





TT 



Table 1. The ranges of the parameters n,//, Sq. 



P := 1/T ^ 1^1, and we take the T = 0-case as a point of reference for the ranges 
n, fj,, B in the low-temperature regime. To avoid possible singularities, we assume a 
finite non-zero density < n < oo, i.e. — 1 < ^ < 0, < iJ < tt. Thus the two zeroes 
of e are supposed to lie inside the interval tt, 7r[. 

Since £(|A:| > i?) > 0, terms ^ exp \— (3£(±tt)] are exponentially small. This allows 
us to eliminate exponentially small quantities by performing two integrations by parts 
in Eq. (6): 

^-rin(l + e^-W) A'^ 



" £(fc')fc'lfe -/fee(fc')dfc' ,,,,dfc 1 e(fc)fc-/^e(fc')dfc' 



4cosh^^^ 27r 27r i + e-^'^('=) 



2 



(8a) 
(86) 



fc+(«) rk+{u) 



7-00 4cosh2| 2tt' 

As explained above, the boundary terms stemming from the integrations by parts 
are exponentially small and are therefore neglected. In the last line, u = 0e{k) was 
substituted and the integration boundaries sent to ±oo, which results again only into 
exponentially small corrections. Since u/P is small by construction, we Taylor-expand 
s{k) around ±B to determine the boundaries k± for the integrand in Eq. (8c): 

^ = {k+- B)e\B) + i(fc+ - Bfe"{B) + \{k+ - Bfe^^\B) + 0{kX) 

_ u e"{B) 2,{e"{B)Y-e'{B)e^^){B) u^ , 

+ ~ ^pe'{B) 2£'3(B)/32+ 6£'5(S) /J3 + ^ l^^ )■ 

Analogously, 

. _ „ u e"{B) 3{e"{B)r-s'e(^HB) u^ 

(3e'{B) 2e'3(B) /32 Qe'^{B) ^'^ >' 

We now expand the integrand of the outer integral in Eq. (8c) in powers of u. Only 
even powers contribute, leading to integrals of the type 

1 



/. 



oo 



4 cosh' 



2 



■ dw = 1 (9a) 



2 



oo 2n 

du = (2n - l)!4n (l - 2^-^^) C(2n) , n = 1, 2, . . . (96) 

4i 0L)i3li 2 
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Especially, C(2) = tt'^/G, C(4) = Tr^/gO. Then 

In order to identify the T-dcpendence of e{k) and B, one performs an asymptotic 
low-temperature expansion of the convolution in Eq. (7) in an analogous way: 

J —: 



K{k -k')\n(l + e-^'^'^/^) ^ = -/? / K{k- k')e{k') 



27r 



+T 



12 £'(B) 

Thus the integral equation for e{k) takes the following form: 

e{k)=e^°\k)+ n{k - k')e{k') — 
J-B 27r 

where (we drop the neglected O (T^) in the following) 

£(0) (fc) = _ cos A; - + K{k- k')e{k') — - T 
J-B 27r 



(11) 



TT k(A: - B) + k(/c + S) 



_B z/i 12 e'(B) 



-i!^ / rial"! - jsai -m^Mt^m-m 



Consider now the function p(fc), defined by 

^K(fc-fc')p(fc')^ (13) 

p(°)(fc) = ^. (14) 

Note that this function is equal to the density of Bethe ansatz roots in the 
thermodynamic limit only for T = 0. For T > 0, the integration boundary B acquires 

a T-depcndence, and so docs p. 

Due to the symmetry K{k) = K{—k), the equality J^^ p^°\k)s{k)dk = 

J-B p{k)s^'^^ (fc)dfc holds. Using this relation and Eq. (12), we can manipulate Eq. (10) 
further and obtain 

.2^' PiB) 



g{T,p)= [ p{k){-cos{k)-n)Ak-T%P] 

J-B -J ^ 



-B 

15 



(B) 



l(e"(S))2 e(3)(B) \ p'{B)e"{B) p"{B) 



(15) 



^4 e'^'{B) 12e'4(5)y Ae'^{B) 12e'^{B)_ 

We now have to identiiy the T-dependence of p{k), s{k) and B. Since in the equations 
for p{k) and s{k), only even powers of T occur, we make the T-dependence explicit as 



B = 
p{k) = 

e{k) = 
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Bo + BiT^ + B2T^ (16a) 

poik) + piik)T^ + P2ik)T^ (166) 

£o(fc)+ei(fc)r2 + £2(fc)T^ (16c) 

Let us consider Eq. (11) first. After inserting (16a, 16c), one obtains equations that 
determine e^, u = 0,1: 

dfc' 

f-Bn 



£o(A:) = — cos k — p ■ 



Bo 



K{k- k')so{k')- 



TT r ° dk' 

eiik) = --—[K{k-Bo) + Kik + Bo)]+ K{k-k')ei {k') — 

^^^0 J-Ba 

where we abbreviated eg(i?o) =: £o- Obviously, B\ is given through 
exiBo) 



(17a) 
(176) 



Bi = 



£[,(Bo) 



(18) 



Note that _Bi > 0. It will become clear below that we do not need ei{k). Analogously, 
equations for px^i are derived: 



/Ba J I./ 

«(A:-fc')p>'(fc')^ 
-Bo 27r 



(fc) = ^ \K(k - Bo) + K{k + Bo)] 



Po 



-2n 



P^^\k) = ( B2P0 + B^px + ^B'ip'^ ] [K{k - Bo) + K{k + Bo 



2-K 

+ ^BIpo [K'{k - Bo) - K'{k + Bo)] 
47r 



(19a) 
(196) 

I 

(19c) 



where p„{Bo) =: Pu etc. (see Eq. (14) for po)- 

On the other hand, with ds{k) := — cos k — p, 



i: 



p{k) [— cos k — p]dk = 



po{k)ds{k)dk 



pBo 

/ Pi{k)d, 

J -Bo 



{k)dk + 2BiPode{Bo) 



(20) 



P2{k)de{k)dk + 2 



B2P0 + Bipi + -Bfp'o 



d,{Bo)+BfpodUBo) 



We insert now (2.2) into Eq. (20) and make use of eo(-Bo) = 0. One ends up with 

/B /"Bo 
p{k) [- COS k - ij]dk= / po{k)[-cosk- p]dk + T^poeoBf. (21) 
-B J-Bq 



Let us look at the second term on the rhs of Eq. (15). Here, we identify a T^- 
contribution as well, originating in: 



em e'o ^ 



./2 



fp'. 


4pq\ 


\e'o 


e'i )\ 



(22) 



where, again, £j/(i?o) =: £1/ etc. In the last equation, it is convenient to represent p\ 
in terms of ei, which follows from Eqs. (176, 19a, 196): 



pi{k) = ^Bipoe'oe\{k). 



(23) 
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Taking Eqs. (15, 21, 22, 23) together, we arrive at the low-temperature expansion of 
9- 

g{T,f,) = 5o(m) +51(^)^2 +52(m)T^ + O (T^) 

dfc 



50 (m) = 



51 (m) 



Po{k) (— cosfc — n) 



52(m) = - S-rr 




277 



(24a) 
(246) 



(24c) 



if ^ 
124''^ 



+ 



I2e'i 



(24d) 



The order has been obtained previously [11] and is common to one-dimensional 
critical systems [20]: The velocity Vc associated with the low- lying excitations reads 
Vc = £q/{2t:po), so that gi{p) = Let us shortly comment on g^- In the first lino of 
Eq. (24d), only eQ,po and derivatives of these functions taken at Bq enter. All these 
quantities are obtained from the T = Bethe ansatz, analogously to v^- Additionally, 
the second line of Eq. (24(f) contains ei, given through Eq. (176). 

From (2.2), the charge susceptibility Xc{T,iJ,) = —d'j^g{T, jj) as a function of T 
and jji is derived. When calculating the specific heat however, attention has to be 
paid that it is defined at constant particle density, whereas above, calculations have 
been carried out at constant chemical potential. Thus the thermodynamic potential 
acquires a T-dependence [21] 



dxPln = 

and we have to consider 
C{T,n) 



dij.n\T 



du,n\T 



(25) 



(26) 



Analogously to (16a, 166, 16c), we write the T-dependence of fj,{T,n) explicitly, 

M(T,n) =/xo(n)+/xi(n)T2. 
Then, according to Eq. (25), 

1 7r„ 1 



^i(n) 



-da 



6 



(27) 



6xc(T = 0,n) 

where Xc{T = 0,n) = Kc/vc, with Kc being the Luttinger parameter, to be obtained 
from Kc = 4:Trpo{Bo) [11]. Note that the same symbols Xc, Vc are being used for the 
susceptibility and velocity as functions of /ip and of n. The T^-dcpcndcncc of p leads 
to an additional contribution to the T^-term in C/T, stemming from the first term 
on the rhs of Eq. (26). The final result is 



C(T,n) 



n 



-12sf2(Ato(n)) + 



da 



Xc(T = 0,n) 



(28) 
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Let us shortly comment on the integrable Lieb-Liniger Bose gas, Eq. (5). The 
analogues of Eqs. (6, 7, 13) read 

1 fB Q 17/ 

,<,(r,ri.-T£,.,(i + .-.*.)f , 

Setting up equations for po,i, £o,i as above, the results (246-24d) yield 5o,i,2- The 
case of very strong interaction c/n 1 allows for the explicit calculation of go, 1,2 

inchiding the order n/c, because in this case, the integration kernel in Eq. (29) (and 
in the equations for eo,i) becomes independent of the spectral parameter. As a result, 

1 / 2Bo\ „ 2n^7r 
Po{k) = — 1 + — -], Bo = mr+ 



27r \ CTT 

e'oik)=2k, .r(fc) = -^ 
2n^Tr^ / 6n7r\ 

_ 1 1 

~ ~ 12n ~ 3c 

7 11 

92 



3c 



where n = n{p) with ^{n) = n^ir"^ — ^^^J' ■ The chemical potential and the specific 
heat at constant density read 



C(T,n) /I 2 \ / 1 



6n 3c y Vl2n57r2 AOn'^n^c 



2.3. The q = 0-case 



In this section, we focus on the case g = in Eq. (2), i.e. 77 — > 00. The corresponding 
model has been termed the "phase-model" in [15, 16]. This case is intriguing in 
so far as ground-state properties at T = can be calculated explicitly (i.e., the 
corresponding linear integral equations can be solved), while at the same time, the 
model still describes correlated bosonic particles. The reason for the solvability of the 
T = 0-cquations is that the integration kernel K{k) = 1 here. Since in the previous 
section, the coefficients of the low-temperature expansion have been expressed in terms 
of T = 0-quantities only, these coefficients can be calculated explicitly for g = 0. 
Furthermore, a very accurate numerical iteration procedure is possible in this case, 
which allows us to test the analytic results against the numerics. 

We defer the detailed calculations to the appendix. The specific heat at constant 
density reads: 

C{T,n) ^TT n+l ^3 (n+l)(23 + 9cos|^) 
T 3 sin ^ + 60 sin^^ 
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^5 (n + 1) (33183 + 28012 cos |^ + 1525 cos 



+— — 5 — 'lajLT^ + o(T^). (30) 

12096 sin^^ -r \ ) \ ) 

It is worthwhile noting that all coefficients are positive. Furthermore, for n — > 0, the 

coeiEcient of the order T^'' diverges ^ n~*'^+'^, whereas for n oo, the divergence 
is ~ n^''+^. Note that both divergencies are excluded from the beginning by having 
restricted to — 1 < /it < 0, that is, the density is supposed to be finite non-zero. 

3. Comment on the high-temperature expansion of the g-Bose model 

In this section, we indicate how to perform a high-temperature expansion of the q- 
Bose model at fixed density n. The general formulation is presented first, before 
concentrating on the q = 0-case. The goal is an asymptotic expansion of the form 
g{T,fj,)/T = gu{lj)P^ (in a sUght abuse of notation, we use the same symbols for 
the coefficients as in (1)), with = 0, 1, 2, 

3.1. General formalism 

Let us write the integral equation Eq. (7) in terms of Ina(fc) := —(3s{k) and A = 1 + a: 

dk' 

K{k-k')\nA{k') — . (31) 
27r 

The /c-dependent driving term on the rhs is bounded. Similarly to the low-temperature 
expansion, we make the T-dependence explicit in the following ansatz: 

lnA = ao + ^ + ^ + 0{T-^) (32a) 

Ina = In (e- - 1) + T'^ + '"' + -"^^7^ + '""^ T"^. (326) 

^ ' l-e-«o 4cosh2^^ ^ ' 

Since n is to be held constant, we have to retain /x in Eq. (31) so that its T-dependence 
can be determined. Thus in zeroth order, 

ao = -ln(l-e'^^) . 

Furthermore, from Eq. (6) it follows for the particle density n that n = d^^aQ, which 
results in 

71 

/3M=ln— -. (33) 
n+1 

Here, the leading T-depcndence of such that n is constant is displayed : /i ~ T. 
From this, the canonical potential gc '■= g{T, i-i{T, n)) + ii{T, n)n is calculated: 

gc= -T[{n+l)\n{n+l) -n\nn] (34) 

The corresponding entropy 

S{T, n) = (n -I- 1) ln(n -|- 1) - n In n (35) 

can be understood as follows: Consider N (classical) particles to be distributed on L 
lattice sites. There are 

possibilities of doing this. The associated combinatorial entropy per lattice site is 
Sq{N,L) = j^\nil{N, L). In the thermodynamic limit N,L ^ oo, N/L = n fixed, 
this agrees with (35). 
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Proceeding further, one finds the linear integral equation for the coefBcient a^: 

c\k' 

K{k-k')ai{k') — , (36) 
-TT 27r 

and the corresponding expansion of g: 

dk' 

ai{k') — . 

From Eq. (36) it follows that J^^ ai{k')dk' = 0, such that the leading T-dependent 
contribution to is O (/?) . 

In order the insertion of Eq. (3.1) into Eq. (31) yields 

-alik) + 2a2{k) (e"» - 1) 



4cosh''if 



n{k-k')a^{k')^. 



Once a2 is determined from this equation, one has 



^=ln(l-e-'^'^^ 



a2{k')^ + 0{r') 



from which one can calculate n = —dfj^g{T, /i) and therefrom the T-dcpcndcncc of /i 
for fixed n, i.e. the coefficient in the expansion ii{T) = /kq + (3 Hi- Eventually, one 
arrives at the canonical potential gc{T, n) := g{T, ij,{T, n)) — iJ,{T, n)n, from which one 
obtains C(T)|„ with the leading T-dcpcndcncc C(T) ^ /?^. We do not carry out this 
program here in its full generality, but rather concentrate on the q = 0-case, where 
the integral equations can be solved and thus the coefficients of the high-temperature 
expansion calculated explicitly. 



3.2. The q = 0-case 

For g = 0, we take advantage of Eq. (A.l) in Appendix A. By expanding the integrand 
there in powers of (3cosk, one obtains 



ln(l + e") + 



13' 



P"^ -2 + cosh a; 



16cosh^f 512 cosh^f 



33 - 26 cosh(a;) + cosh 2x 



73728 



cosh 



6 X 



J x=/3(/i-s(T,^)) 



Furthermore, also from Eq. (A.l), the particle density is obtained as 



a(T,M)=£- 



dk 



^g-/3cosfe-/3/i+/3g(T,^) 2^' 

which yields n in terms of /i — g{T, fi): 



(3- 



X ^ 3 + cosh x + 2 sinh x , x 

e tanh - ^ e ^ tanh - 

4 2 128 coshes 2 



(3' 



18432 



42+16cosha;-|-10cosh(2a;) + 54sinha; + 9sinh(2a;) , x 

X ^—^ i — - tanh - 

cosh f 



J x=0{iJ.-g{T,ij.)) 



(37) 



(38) 
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Figure 1. The specific heat for q = 0, n = 0.5, 1, 2, 4, together with the low- and 
high-temperature asymptotics according to Eqs. (30, 39). 



From Eq. (38), — g) is obtained as a function of n order by order in This 
result is then used in Eq. (37) to express g in terms of n. Finally, the first terms of 
high-temperature-expansions of the canonical potential g^ and the specific heat read 



gc{T,n) ^ -T[{n + l)\n{n + 1) - nlnn] - /3 



4 n{l + n^) 



4(1 -hn) 64(1 -hn3) 



6n 



2^4) 



C(r,n) = 



2(1 



1152(1 + n)5 

_ 3n(l + n^) ^, ^ 5n(2_ 



16(1 



192(1 -hn)5 



(39) 



5.5. Comparison with numerics for q = Q 

The difficulty in a numerical calculation lies in the steep (~ (3) gradient of e{k) around 
k ~ ±B. The induced numerical inaccuracy generally is too high to allow for a 
quantitative check of the coefficients (although the exponents in the low-T-expansion 
could be identified, [11]). However, at q = the iterative numerical approach becomes 
particularly accurate because only one non-linear integral equation has to be solved 
(see Eq. (A.l) in the appendix). Results for different n at g = are shown in figures 
1, 2. 



4. Low-temperature asymptotics of the XX Z-raode\ in a field 



It turns out that a low-temperature expansion for the free energy of the spin-1/2 
XXZ-m.ode\ in a magnetic field is technically very similar to the one given above for 




Figure 2. Low- and high-temperature asymptotics of the specific heat at q = 
and n = 0.5, 1, 2, 4. In a), the order of C(T) is shown on a double logarithmic 
scale, after having subtracted the T, T^, -contributions according to Eq. (30) 
(denoted by C^^^lt in the figure). The dashed straight line is oc and a guideline 
to the eye. Analogously, in b), the /^^-contribution in C is shown, where the 
orders /3^,/3*,/3^ have been subtracted following Eq. (39) (denoted by Csx,ht in 
the figure). Here, the straight dashed line is oc 



the g-Bose model. This is not too surprising, as argued in the introduction. 
The Hamihonian is defined as 

L L 

H = jJ2 [SiS^+, + SpJ^, + (cos 7)^1 ^l+i] +hY,Sl, (40) 

1=1 1=1 

where < 7 < tt is chosen such that the model is critical and low-energy excitations 
have a linear dispersion relation. The summation runs over the L lattice sites of the 
chain, where periodic boundary conditions are imposed. 

Currently, three (equivalent, [22]) approaches to calculate the thermodynamics 
exactly are known: the TBA, the QTM and the single-integral-equation approach 
[23]. For our purposes, the formulation developed within the framework of the QTM- 
technique is best suited. The free energy per lattice site is given by 

/(T,/i) = eo-|^- [d*ln»S] (0) (41) 

ZTT 7 

lnb(a;) = - l3vfU{x) + — + [k * In «8 - k_ ln»] (x) 

2(7r-7) 

lnb(a;) = - ^vfUix) - — + [k * In <8 - k+ In »] (x) 

2 7r-7 
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,(0) ^ 

*B =l + b, QS = l + b 
d{x) ^ 



cosh -X 

7 



k{x) 



sinh(f -7) fce'*^^ dk 



2cosh^fcsinh(f 27r 

K±{x) = k(x ± i7 =F ie) 

and convolutions are defined as [/ * g] (x) := ,f(x — y)g{y)dy. The quantity eo is 

the ground-state energy, vi^'' the spin velocity at zero magnetic field. We now focus 
on the low- temperature regime, that is, jSh ^ 1. In this case, ln*B ~ exp[— /J/i] 
is exponentially small and thus is neglected. It is convenient to introduce l3s{x) := 
Inb(x). Then 

/(T, h)-eo = -— -d{x) In ( 1 + e'^^^^M da; 

= -TT / -d{x)e[x)dx - — / -d{x) In ( 1 + e-'^^^^M da; (42) 
27r J_oo 7 27r J.^^ 7 V / 

e{x) = -v^°U{x) + -h + T K{x-y) In ( 1 + e^^^A dy 

2(7r-7) J_oo ^ ^ 

2{tt - 7) 

+ / k{x- y)e{y)dy + T I k{x - y) In (l + e-'^^(^)) dy. 
J— 00 J —00 ^ ' 

Using Fourier transform techniques, the last equation is written in a more compact 

manner as (we rename the spectral parameter as k) 

£{k) = -vf^^kx{k)+l3h+ r k^{k-k')\n(l + e-^'^'''A— (43) 

_ 2sin(27) 

cosh(2fc)-cos(27)' 

which makes the similarity with Eq. (7) apparent. 

One now goes through the same steps as in the analysis above. We give the 
corresponding results in the following. The functions po,i satisfy the linear integral 
equations 



/ dk 
k^{k-k')p.{k')^ 
-Bo 27r 



p(°){k) = l-d{k) 



27r 

and for the functions So,i{k), we find 



p^^\k) = Bi^ [k^{k - Bo) + k^{k + Bo)] 

ZTT 



/Bo Ay 
k^{k - k')e,{k')^ (44a) 
-Bn 27r 
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;(0)(yfc) = _1;(0)2kt(A;)+/i 
TT ^ 

= --^ [k^{k - Bo) + k^{k + Bo)] + [ k^{k - k')ei{k') 

^■^^0 J-Bo 



2^' 



(446) 
(44c) 



The shift Bi is given by Eq. (18), with sifi specified in Eqs. (4). The low-temperature 
expansion of the free energy reads 



/(T, h) = fo{h) + h{h)T^ + h{h)T^ + O (T^) 
k{h) 



= eo + 



d{k)vf''^kx{k)dk 



|fe|>So 



h{h) 

f2{h) 



3 s'o 

4 3 



(45 a) 
(456) 

(45 c) 



Ml! 



,(3) 



1244 



-/3 



(45 d) 



^0 ^0 

where £0 := eo(-So) etc. 

In this model, the particle density is constant per construction, so that the specific 
heat can be calculated directly from C{T,h)/T = —drfiT^h). However, formula 
(26), with replaced by h, can be used to calculate the specific heat at constant 
magnetization m. The results in Eqs. (45an45cf) are confirmed for free fermions, cf. 
Appendix B. 



5. Summary and outlook 

The g-Bose model and the critical spin-1/2 Heisenberg chain in an external magnetic 
field were studied in the regime of low temperatures. In both cases, the thermodynamic 
potentials could be expanded in terms of even powers of T. We developed a method 
that can be used to calculate the coefficients of this expansion, and carried out the 
calculation for the T'*-contribution. 

Directions for future research are twofold: One question is in how far the results 
obtained here can be used to determine the coefficients of the leading irrelevant 
opc;rators in the low-energy cfFcctivc; Hamiltonians for these models. The other 
question concerns the relationship between the spin chain and the bosonic model. 
The absence of backscattering in both models (at low energies for the spin chain) 
suggests that a relationship exists on a more fundamental level. 
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Appendix A. Specific heat at q = 

By combining Eqs. (6, 7), it is not difficult to show that for q = 0, 
g{T,ti) = -T In (l + el3'^osk+0^.-Pa(T,f.)\^ g_ 

In the low- temperature regime /3 » 

(— fccos fc + sinfc) sin /c dk 



^^^'^^ ^ '-^4cosh^{cosk + fi-g{T,n))2n 

which is obtained by two integrations by parts, neglecting exponentially small 
quantities. One now substitutes u = /3(cos k + ^ — g{T, /i)) and obtains 



1 - - (/^ - fl)) - - (m - g)) arccos ^-{^l-g) 



4cosh2f 27r' 



where the integration boundaries have been sent to ±oo, again neglecting exponentially 
small corrections, and the (T, /i)-dependence of g has not been displayed explicitly on 
the rhs. One now expands the integrand in powers of u/ (3 and makes use of Eqs. (9a, 
96), resulting in 

g{T, m) = - - fVl - + (m - 9) arccos(5 - (A.2) 



6 VI - (m - 9? 
JTT^ l + 2{fi-gf 
360 (1 - (a* - ^ ■ ' 

_ 6 3l7r^3 + 8(M-fl)^(3+(M-g)^) . 8^ 
5040 (1 - (m - 9)2)9/2 ^ ^ 

This equation for (/ can be solved order by order by expanding g = go + giT^ + 52?"* + 
gsT^. This allows for the calculation of n{iJ,,T), and therefrom, by holding n fixed, 
the (T, n)-dependence of /U is derived: 

/z(T, n) = - — + - sm Bq 

\ n + 1 TT 

^ 7r(-l-n + 7rcotBo) ^2 
6(n + 1) sin So 

TT^ (2937rcosBo + 277rcos(3Bo) - (n + 1) (37sinBo + 9sin(3So)))y4 ,^ 

1440(n + 1) sin^ Bq ^ ' ' 

(10888827rcosSo + 2758857rcos(3Bo) + 111537rcos(5Bo) 



725760 sini° Bq 

-{n + l)(38354sinSo + 26487 sin(3So) + 1525sin(5So))) 

Bo = (A.5) 

n + 1 

Taking this T-dcpcndcnce of fi into account, one arrives at the specific heat, taken 
at constant n, Eq. (30). Furthermore, the charge susceptibility is obtained as 
XciT,n) = l/dnf^iT,n). 
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Let us now confirm these results from the general formulae (24oh24d, 28). 
Eqs. (19a) (for v = there) and (17a) are solved by 

tiiuBii + BoH 
eo(fe) = — cosfe — /X — . 

IT - Bq 

The condition So{±Bo) = yields 

sinSo , ( Bo \ 
IJ, = h 1 cos i?o 

TT \ TT / 

According to Eqs. (7, A. 2), e(fc) = — cosfc — ^ + g, and thus cosBq = go — 

Bo 

-B, 



On the other hand, we know [11] that J Padk/2 TT = n, which yields Eq. (A. 5). 
This confirms the T^-order of Eq. (A. 4) and is consistent with n(/zo) = — 9^o(7o(/^o)- 
However for the sake of brevity, we do not display the /^o-dependence of n. The 
coefficients (24c, 24d) can thus be calculated and read 

TT 

91= - 



92 = 



6wc 

TT f n TT^cosSo 77r3(l + 2cos^ So) 



TT-Bo\ 72t;2sinSo 36vcSm^ Bq 360sin^Bo 



1 . TTn 
Vr = sm ■ 



n+1 n+1 

These results are confirmed by the slightly different analytical calculation sketched in 
the first part of this appendix. Note that the q = 0-case is special in so far as po {k) is 
a constant, as well as £i{k). 

Appendix B. Specific heat for free fermions 

We calculate the first terms of a low-temperature expansion of the free energy for 
free spinless fermions on a lattice in two different ways: Directly and using the exact 
solution. Both are equivalent and yield the same results. 

At 7 = tt/2, a Jordan- Wigner transformation followed by a Fourier transformation 
of the Hamiltonian (40) yields 

H = ^{J cos kj + /i)4 . Ckj with kj = —j. 

Thus in the thermodynamic limit, 

-/3/(T, h) = — J \n[l + e-'5('''^°«^+^) j dx, 
which is very similar to Eq. (A.l). Proceeding analogously to there, one finds 

/(T,/i)= I Wl_ ( :^ ) -l^arccos- 



TT 




J J 



-J 



iTT^ (TV l + 2{h/jf 



(B.l) 



360 V J / /T (un\r >/'^ ' 
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On the other hand, we may use the exact solution Eq. (4) with 

p(k) = 

IT cosh 2x 

1 .J 

B = -arccosh— . 

2 h 

Note that p{k) = po{k), s{k) = So{k), B = Bq here. Inserting these values into 
Eq. (4), one confirms the result (B.l). 
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